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1. Introduction 



A puzzle of continuing interest in particle physics is the apparent quantization of 
electric charge. The electric charge of all the particles we see appears [1] to be an 
integer multiple of Qe/S, where Qe is the electron charge. However it is theoretically 
consistent to have particles with electric charge eQe where e is any real number. We 
consider e < 1 in this paper, and refer to these particles as "milli-charged." 

Milli-chargcd particles can be introduced into the Standard Model in a variety of 
ways. One is to add an SU{?)) x SU (2) singlet Dirac fermion with hypercharge Y = 2e. 
However, this is not simple [2] if the hypercharge U{1)y is embedded in a grand 
unified gauge group. A second way of generating effectively milli-charged particles, 
which works even if hypercharge is quantised, is due to Holdom [3] . He introduced a 
second unbroken "mirror" C/(l)', and showed that the photon and the "paraphoton" 
can mix, so that a particle charged under the 11(1)' appears to have a small coupling 
to the photon. A third mechanism for introducing milli-charged particles is to remain 
with the Standard Model particle content and allow the neutrinos to have small 
electric charges [4]. If the Standard Model hypercharge operator Ysm is redefined to 
be y = Ysm + 2 ei(i?/3 — Lj), neutrinos acquire small electric charges and the 
Standard Model anomaly cancellation is preserved. We assume here that U{l)em. is 
an unbroken symmetry; it can be spontaneously broken in models with milli-charged 
scalars [20], but the bounds in this case are slightly different from what we discuss 
here [21] (for instance the photon has a mass). 

The purpose of this paper is to collect and update experimental and observa- 
tional bounds on milli-charged particles. Such limits have been previously considered 
for milli-charged fermions that arc not neutrinos [6, 7, 8, 9, 10, 11, 12, 13, 14] and for 
milli-charged neutrinos [5, 15, 16, 17, 18]. See also [19] for a discussion of astrophys- 
ical constraints on all types of milli-charged particles. We recalculate the laboratory 
bounds using recent data, and include the limits from recent experiments searching 
for miUi-charged particles [27, 28]. We numerically evaluate the constraint from Big 
Bang Nucleosynthesis on particles with milli-hypercharge, and find a bound similar 
to previous estimates. We also revisit the astrophysical constraints, resolving a dis- 
crepancy in the White Dwarf bounds between [9] and [12] (in favour of the stronger 
bound of [9]). Our Supernova limit differs slightly from what is found in Rcf. [10], as 
will be explained later, and the Red Giant bound remains unchanged. We quote the 
previous estimated bounds [12] from balloon experiments searching for strongly in- 
teracting Dark Matter, and from underground WIMP detectors. For a more detailed 
discussion of the calculations, see for instance [12, 14, 17, 19]. 

Our conclusions are presented in Fig. 1 which is an update of Fig. 1 of Ref. [14]. 
It presents bounds on milli-charged fermions in models with an extra U{iy, and 
without. MiUi-charged neutrinos will be briefly discussed at the end of the paper. 
The hmits from laboratory experiments, from Red Giants and White Dwarfs are 
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Figure 1: Regions of mass-charge space ruled out for milli-charged particles. The solid 
and dashed lines apply to the model with a paraphoton; solid and dotted lines apply 
in the absence of a paraphoton. The bounds arise from the following constraints: AC — 
accelerator experiments; Op — the Tokyo search for the invisible decay of ortho-positronium 
[27]; SLAC— the SLAC milli-charged particle search [28]; L— the Lamb shift; BBN— 
nucleosynthesis; CI — < 1; RG — plasmon decay in red giants; WD — plasmon decay in 
white dwarfs; DM — dark matter searches; SN — Supernova 1987A. 



independent of whether there is a paraphoton, and appear as solid lines. We discuss 
the Supernova and nucleosynthesis bounds with some care for the model without 
a paraphoton; the limits should be similar in the presence of a paraphoton [14] so 
these bounds also appear as sohd lines. The hmit from < 1 diflFers in the two 
models: in the absence of a paraphoton, the millicharge annihilation cross-section is 
proportional to e^/m^. So the excluded area, indicated with dotted lines, extends to 
lower masses for smaller e. The bound on models with a paraphoton is the horizontal 
dashed hne. 
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2. Laboratory Bounds 



We plot four lines on Figure 1 from laboratory data: 1) a combined "accelerator" 
line consisting of the limits from LEP and beam dump experiments [7, 12]^, 2) the 
bound found by the Tokyo group [27] from the non-observation of invisible Ortho- 
positronium decay, 3) the limit from the dedicated milli-charged particle search ex- 
periment at SLAC [28], and 4) an updated constraint from the Lamb shift using 
more recent data [24]. 

LEP has taken many years of data since limits from LEP were previously consid- 
ered [12], so we briefly discuss possible bounds, although these are weak because the 
milli-charge coupling to the Z is suppressed by sin^^^K- A search for particles with 
fractional charge e = 2/3 was performed by OPAL using 1991-93 data [22], which 
rules out e > 2/3 for < 84 GeV. This bound could be extended to the present 
kinematic limit < 100 GeV, if one assumes that a particle with 1 > e > 2/3 would 
be seen as such in the detector. Fractionally charged particles would contribute to 
the invisible width of the Z if they were not seen in the detector ^, in which case the 
LEP bound can be extended to 

e < 0.24 m, > 45 GeV, (2.1) 

from requiring that that milli-charges not contribute more than the 2(7 error to the 
invisible width of the Z at LEPl. 

We calculate an improved bound from the Lamb shift measurements by requiring 
that the milli-charged particle vacuum polarisation contribution be less than the 
2(j error ( = 20 xlO~'^ MHz) [24]. This is a more recent measurements of the 
2-S'i/2 — 2Pi/2 Lamb shift [24] than used in [12]. There are more precise measurements 
of combinations of Lamb shifts of different principle quantum number n {n ^ 2) [25]; 
however they do not substantially improve our bound so we use the 2Si/2 — 2Pi/2 
results because the milli-charged contribution is easier to disentangle. The bound 
from the Lamb shift is stronger than from g — 2, despite the fact that — 2 is measured 
very precisely [26], because milli-charged particles contribute to the Lamb shift at 
one loop but to (jf — 2 at two loop. For a discussion of the effects of milli-charged 
particles on precision QED measurements, see [9, 23]. 

3. Big-Bang Nucleosynthesis 
3.1 Dominant Processes 

Particles with small electric charge will interact with the plasma in the early universe 

^The values of excluded by direct searches are quoted in Table II of [12], but the values given 
are in MeV, not GeV as claimed in the caption. The limit is placed correctly, however, in Fig. 2 of 
[12]. 

^This bound assumes that particles with 1/4 < e < 2/3 would look like noise in the detector, 
and not be mis-identified as e = 1 tracks. 
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and in this way get thermally excited. Therefore, their primordial energy density and 
thus their interaction strength can be constrained by the standard nucleosynthesis 
hmit on A^eflf, the effective number of thermally excited neutrino degrees of freedom. 
A milli-charged fermion must be a Dirac particle and thus has four intrinsic de- 
grees of freedom. If it were fully excited in the early universe, it would contribute 
AA^eff = 2 prior to the annihilation of electrons and positrons ^. Such a particle 
would be excluded by a large margin. During this annihilation process entropy is 
transferred to the milli-charged particles, heating them relative to the neutrinos. 
At the same time, the photon temperature, Ty, relative to the neutrino tempera- 
ture, T,y, is less than in the standard model. Normally, T„/T^ ^ 0.714, whereas 
here T^/T^ ~ 0.849. Altogether, neutrinos plus milli-charged particles contribute 
A^eff = 3 X (0.849/0.714)4 + 2 x 1.401^ = 13.7, so that i^N^^ = 10.7 at late times. 

While it is easy to estimate AAeff for milli-charged particles where full thermal 
equilibrium obtains, we can derive a more restrictive limit by studying its value as 
a function of e in a regime where full equilibrium is not achieved. To this end we 
assume that the new particles interact only electromagnetically by their milli-charge. 
The bounds in the case where there is an extra U{1) are at least as restrictive, due 
to the presence of the paraphoton [14]. 

We find that the dominant electromagnetic processes which couple the gas of 
milli-charged particles / to the electromagnetic plasma are 

7 - ff, 

ef ^ ef. (3.1) 

The rates for these reactions are proportional to while those for other conceivable 
processes such as 77 ff or 7/ — 7/ are proportional to and thus much 
smaller. At low temperatures, the electron-positron number density is exponentially 
suppressed, making the // production rate approach zero quickly when T ^ mg. 
The third process in Eq. (3.1) does not contribute to the production of milli-charged 
particles but helps to achieve kinetic equilibrium. 

For a calculation of the // production rates we use Boltzmann statistics for all 
species. Then the rate for the pair process is [38] rg+g-_,jj = ne{cr\v\), where {(j\v\) 
is the thermally averaged annihilation cross section. It can be expressed as [36] 

1 f°° 

(''W) = SruiTKUruJT) L. * ^" " acM, (3.2) 

where is a modified Bessel function of the second kind and iTcm is the center-of- 
mass cross section. The squared matrix element is 

- J2 = _ [iPe-Pf){Pe+Pf) + {Pe-Pf-){Pe+Pf) + m){pe-Pe+) 

^This is when the BBN bound on the number of neutrinos applies 
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+ ml{pfPf) + 2mlmj , (3.3) 

which translates into the usual CM cross section ctcm = (47r/3) e^a^ if all the 
particles can be considered massless. In this limit we find 

(re^e-^ff) = ^ eWT, (3.4) 

where C refers to the Riemann zeta function. 

The decay rate of a transverse plasmon (photon) with energy cj and momentum 
k into massless milli-charged particles is ^j^ff = e^aZiuP' — k'^)/'iuj [19]. In a 
relativistic plasma, the plasma frequency is given by a;| = (47r/9)Q;T^. Except for 
the low-energy part of the blackbody photon spectrum we have cj ^ cjp, a limit where 
the dispersion relation is u'^ — k'^ = |c(jp and the wave-function renormalization factor 
is Z = 1 [19]. With Boltzmann statistics we find (a;~^) = (2T)~^ so that finally the 
average plasmon decay rate is 

(r.-.//-) = I eWT. (3.5) 

Comparing the rates for the two processes we find 

implying that the plasma process is actually somewhat more important than pair 
annihilation. Notice, however, that the plasma process is only important if m ^ 
ujp/2. During BBN, the plasma frequency is roughly given by cjp ~ ^J (47r/9)aT ~ 
0.1 T, so that at T = 1 MeV, where BBN begins, the plasma process is only important 
for m ^ 50 keV. Therefore the plasma process is important only for milli-charged 
particles with masses in the keV region, where astrophysical bounds from stellar 
evolution are much more stringent. Thus, it seems safe no neglect the plasma process 
in the BBN calculations. 

3.2 Solving the Boltzmann Equation 

The standard procedure for calculating the evolution of number and energy density 
of a given species is to use the Boltzmann collision equation. It describes the time 
evolution of the single particle distribution function, /i, of any given species. In the 
expanding universe, formally this equation can be written in our case as 



^-i/p^ = (ann + Cel)[/l], (3.7) 

where H = R/R is the Hubble expansion parameter. On the right-hand side, Cann is 
the collision operator describing pair annihilation. Cei is the elastic-scattering term 
which includes the combined effect of scattering on electrons, positrons, and photons. 
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The e+e process is of the form 1 + 2— >3 + 4so that the coUision term can be 
written generically as 



Cann[/l] = ^ / d%d%d%A(n, h, h, h) E W?^\Pr + P2 - PS - P4)(2n)\ 

(3.8) 

Here, (fp = d^p/[{27r)^2E], X^|M|^ is the squared matrix element, and Pi is the 
four-momentum of particle i. The phase-space factor is A = fsfill — /i)(l — /2) — 

/l/2(l-/3)(l-/4). 

The dominant elastic scattering term from scattering on electrons and positrons 
suffers from the usual infrared Coulomb divergence. It can be regulated by including 
Debye screening effects in the medium. However, instead of treating elastic scattering 
in any detail we calculate the evolution of the milli-charged particle ensemble for two 
extreme cases, a) No elastic scattering, and b) Elastic scattering is assumed to be 
efficient enough to bring the milli-charged species into complete kinetic equilibrium 
at all times. It will turn out that the difference between these two cases is quite 
small, justifying our neglect of a detailed treatment of elastic scattering. 

The dynamical evolution of the cosmic scale factor is governed by the Friedmann 
equation and the equation of energy conservation [38], 



= (3.9) 

where p is the energy density and P the pressure, both including the effect of milli- 
charged particles. 

We have solved these equations to find the energy density of milli-charged par- 
ticles. We have always neglected their mass, an approximation which is accurate for 
mg up to the electron mass. For larger masses the milli-charged particles can pair 
annihilate into electron-positron pairs at low temperatures. We have always taken 
a zero initial population of milli-charged particles, an assumption which is valid if 
they interact only via their electric charge. 

Adding energy density to the cosmic plasma around the epoch of e^e~ anni- 
hilation perturbs Big Bang Nucleosynthesis [38]. In order to calculate constraints 
on e we have used the nucleosynthesis code of Kawano [32] , modified to include the 
energy density in milli-charged particles as well as the change in neutrino-to-photon 
temperature ratio. The effect of the milli-charged particles corresponds to additional 
neutrino degrees of freedom of 

AA^ 69 X 10^^ ^ X I ^ elastic scattering, _^ 

I 1.39 full equilibrium, 

an approximation is valid for e ^ 6 x 10~^. This is more than sufficient for our 
purposes since all values of e higher than this produce such a large AAefr that they 
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arc excluded by a huge margin. As expected, the effect of elastic scattering on 
the energy density in milli- charged particles is quite small as this process does not 
produce additional particles. 

3.3 BBN Limits 

There are still some unresolved issues regarding the observationally determined values 
of the primordial light-element abundances. For the past few years there have been 
two favoured solutions, namely the so-called High- Helium/ Low- Deuterium and the 
Low-Helium/High-Deuterium solutions [33]. There seems to be growing consensus 
that the first is the correct one. Nevertheless any bound derived from nucleosynthesis 
should be used with some caution since the question of primordial deuterium and 
hehum abundances is not yet fully settled. 

In the present paper we shall use the data on primordial helium obtained by 
Izotov and Thuan [34] and the deuterium data from Buries and Tytler [35] 

Yp = 0.244 ± 0.002, (3.12) 
D/H = (3.39 ± 0.25) x 10-^ (3.13) 

where errors are estimated la uncertainties. These data give the High-He/Low-D 
solution and are completely consistent with standard BBN for a baryon-to-photon 
ratio of 77 = (5.1 ± 0.3) x 10~^° [35]. Used together, these data provide the constraint 

iVeff = 2.98 ± 0.33 (90% C.L.), (3.14) 

on the effective number of neutrinos, thus tightly constraining any non-standard 
nucleosynthesis scenario. 

Together with the less restrictive case (no elastic scattering) of Eq. (3.11) BBN 
then implies a limit 

e < 2.1 X 10"^ (3.15) 

This bound applies to masses in the regime ^ mg. While our result is very similar 
to what one finds from simple dimensional estimates, our calculation is quantitative 
and the errors are controlled. 

4. Stellar Evolution 
4.1 GlobulEir Clusters 

New low-mass particles will be produced in the hot and dense medium in the interior 
of stars and subsequently escape. This new energy-loss channel leads to observational 
modifications of the standard course of stellar evolution and thus can be used to set 
limits on the particle's interaction strength. For milli-charged particles, limits have 
been set from red giants [5, 9, 12, 14, 13, 39], horizontal- branch (HB) stars [19], white 
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dwarfs [9, 12], and supernova (SN) 1987A [10]. For small masses of the milli-charged 
particles, the most restrictive limits arise from HB stars and low-mass red giants in 
globular clusters. 

At the end of the main-sequence evolution of normal stars, the hydrogen in 
the inner part has been consumed, leaving the star with a core consisting mainly 
of helium. In low-mass stars, this helium core reaches degeneracy before it is hot 
and dense enough to be ignited. This process is very dependent on density and 
temperature, and even minor changes in these quantities produce observable changes 
in the brightness at the tip of the red-giant branch in globular clusters. Therefore, 
the core mass at helium ignition as implied by the color-magnitude diagrams of 
several globular clusters implies a limit on any new energy-loss channel. After helium 
ignition, the stars move to the horizontal branch where they burn helium in their core. 
A new energy-loss mechanism will lead to an accelerated consumption of nuclear fuel, 
shortening the helium-burning lifetime which can be "measured" by number counts 
of HB stars in globular clusters. When applied to a new energy-loss channel, usually 
one of these arguments is more restrictive [19]. For example, a putative neutrino 
dipole moment will add to the efficiency of the plasmon decay process 7 — > z/z/ 
and thus enhance neutrino losses. The helium-ignition argument provides far more 
restrictive limits than the helium-burning lifetime argument because the plasmon 
decay is more effective in the degenerate red-giant core. On the other hand, axion 
losses by the Primakoff process are more effective in the nondegenerate cores of HB 
stars so that the helium-burning lifetime argument yields more restrictive limits. 

The emission of milli-charged particles is a special case in that both arguments 
yield comparable hmits of about [19] 

e < 2 X 10-^^ (4.1) 

The reason is that the rate of the plasma decay process 7 ^ // for milli-charged par- 
ticles is proportional to Up, as opposed to the magnetic-dipolc case (cjp) or standard- 
model neutrino case (cup). The low power of the plasma frequency implies that the 
emission rate per unit mass is almost independent of density. This, in turn, implies 
that the core expansion caused by helium ignition leaves the energy-loss rate per 
unit mass nearly unchanged, while it is "switched off" for the dipole-moment case, 
or "switched on" for the axion case. 

An average value for the plasma frequency in the core of a globular- cluster star 
before helium ignition is ujp ~ 10 kcV (in the center it is about twice that) while it 
is about 2 keV in the core of a HB star. Therefore, the helium-ignition argument 
constrains milli-charged particles with masses up to about 5 keV. 

Of course, this sort of argument applies only if the interaction strength is small 
enough that the miUi-charged particles escape freely once produced in the stellar 
core; one can check that this is the case for e ^ 10~^ [19]. For larger charges, the 
particles would contribute to the transfer of energy rather than carrying away energy 
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directly. In order to avoid observable consequences, the efficiency of energy transfer 
would have to be less than that of photons, i.e. the mean free path would have to 
be very short. It is unhkely that there exists an allowed range of miUi-charges e ^ 1 
where all stars would be left unchanged. However, since laboratory limits and the 
BBN argument exclude values for e above 10~^ anyway, a detailed discussion of the 
trapping limit is not warranted. 

4.2 White Dwarfs 

The observed population of hot young white dwarfs is consistent with cooling by sur- 
face emission of photons and by volume emission of neutrinos produced by plasmon 
decay via Standard Model interactions [29, 19]. Blinnikov and Dunina-Barkovskaya 
[29] set a bound on the neutrino magnetic moment, /j-i, < (where Hb — 

e/2me), by requiring that the additional neutrino emission not cool the white dwarfs 
faster than observed. We can translate this bound into one on milli-charged particles. 

For neutrinos or neutrino-like particles coupling to the photon via a dipole mo- 
ment, the energy loss rate can be written as 

whereas for milli-charged particles, it is given by 

aoce^a^^gi, (4.3) 
47r 

where Q1/Q2 is a factor of order unity [19]. Here, we set Qi = Q2 = 1- One then 
obtains 

where eu = 10^^ e and 1112 = 10^^ fJ-/ fJ-B- In order to calculate an average value of this 
ratio over the entire star, we need to perform an "emissivity-average" of the form 

In order to do these integrals, we assume that the white dwarf is a polytrope of index 
n = 3/2, so that 

where r — and p = pcO'^ [37]. For a white-dwarf mass of 0.7 solar masses, this 
gives 
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If wc then demand that the emission rate due to miUi-charges is not larger than for 
neutrino dipole moments, we find 

ei4 < 17. (4.8) 

This bound is consistent with that found in Ref. [9]. 

One may worry that the way the emissivity-average is performed has conse- 
quences for the bound. One could instead use 

qJ Jdrr^Q, 

This yields €u < 15 so that the way the average is performed matters httle for the 
final result. 

The emissivity-averaged plasma frequency is ((10 keV/cup)^) = 0.16 so that 
(cjp) = 25 keV. Therefore, the white-dwarf bound applies to milli-charged particles 
with rUe ^ 10 keV, similar to the red-giant case. 

4.3 Supernova 1987A 

Finally, the stellar energy-loss argument can be applied to SN 1987A where the 
number of neutrinos detected at Earth agree roughly with theoretical expectations. 
If there are other particles contributing to the cooling of the proto neutron star, this 
will reduce the neutrino fluxes and the duration of the neutrino signal. Therefore, if 
we assume that such hypothetical particles freely stream from the core where they 
are produced, one can put an approximate bound on the allowed loss rate [19] of 

(Q/p) ^ 10^9 g-i g-i^ (4 

to be calculated at average core conditions of about 3 x 10^^ g cm^'^ and 30 McV for 
density and temperature, respectively. 10^^ erg g~^ s~^ is the proto-neutron star's 
average rate of energy loss to neutrinos. 

The main production process for milli-charged particles in a SN core is the plas- 
mon process. In a degenerate, relativistic electron plasma it is orders of magnitude 
larger than e+e~ annihilation. The energy-loss rate is [19] 

Qp-^6v(/.^ + ^)t^Q„ (4.11) 

where Qi is again a factor of order unity. Setting Qi — 1, one finds 

e < 1 X 10"^ (4.12) 

This bound matches the one found by Mohapatra and Rothstein [10]. However, they 
considered nucleon-nucleon bremsstrahlung as a production process and notably an 
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(4.9) 



amplitude where the electromagnetic current is coupled to an intermediate charged 
pion. We believe that a naive perturbative calculation of this process in a nuclear 
medium can be unreliable [19]. However, since the bounds are so similar, a detailed 
study of the nucleon process is not warranted. 

If e is much larger than our limit, the milli-charged particles are trapped in- 
side the proto neutron star, and can only escape via diffusion. The main process 
which keeps milli-charged particles trapped is Coulomb scattering on protons ^. The 
differential scattering cross section on nonrelativistic protons is 

^ = 2eV^!il±^ (4.13) 
dQ |q|4 ^ ^ 

where q, with |qp = 2E^{1 — cos6), is the momentum transfer, E the milli-charged 
particle's energy, and 9 the scattering angle. This cross section is strongly forward 
peaked and thus not a good measure for particle trapping: a particle which is de- 
flected by a small angle continues its way out of the star essentially as if it had not 
been scattered at all. Therefore, we rather consider the usual transport cross section 
which includes an additional weight (1 — cos^). Therefore, 

dor ^ {I + cose) 
ail q^ 

is a more adequate measure for the effectiveness of particle trapping. 

In addition, we need to include proton-proton correlations induced by their mu- 
tual Coulomb repulsion, i.e. we need to include screening effects. This is achieved by 
multiplying the cross section with the static structure function S{q) = q^/(q^ + feg) 
so that finally 



C^C^T.eff _ ^2^,2 (1 + C0S6I) 



For nonrelativistic, nondegenerate protons, the screening scale kg is given by the 
proton Debye scale — — AiraUp/T with Up the proton density. It would have 
been incorrect to use the electron screening scale since the background of degenerate 
electrons is much "stiffer" than the protons, i.e. most of the polarization of the plasma 
by a test charge is due to the protons. With these modifications we find a total cross 
section 



o-T,efr - ^2 

where 



(2 + ;^), f2 + z 
m ' 



(4.16) 



^The scattering on electrons was instead considered in [10], but this is less important than 



protons. 



11 



and rje = fJ^e/T is the electron degeneracy parameter. We have used that Up = Ug — 
/Xg/37r^. The transport mean free path is then found to be 



A, 



'T,efr = '^T,eS'np = €^aT z 



(4.18) 



where we have expressed Up in terms of kj^. Taking T = 30 MeV as a typical value, 
and using 1 for the ^-dependent expression, we find that AT,eff exceeds the proto- 
neutron star radius of about 10 km for e ^ 1 x 10~^. 

When e is larger than this, the particles no longer freely escape. Rather, a 
"photo-sphere" for the milli-charged particles is created. For them to be less effective 
at carrying away energy, their transport cross section must be about as large as that 
for neutrinos or larger. Very crudely, we must compare Eq. (4.16) with the weak- 
interaction transport cross section for neutral-current scattering on nucleons. 



where \Ca\ ~ 1.26/2 for protons and neutrons, while \Cv\ ~ for protons and 1/2 
for neutrons. If we use z — 0.01 near the neutrino sphere, the expression in square 
brackets in Eq. (4.16) is about 4. With this value we find 



Taking an average energy of 20 MeV we conclude that e should exceed approximately 
8 X 10~^ for emission of milli-charged particles to be less important than neutrinos. 

While this derivation is somewhat crude, we conclude that milli-charged particles 
in the range 10~^ ^ e ^ 10~^ are excluded. The plasma frequency in a proto-neutron 
star is roughly 10 MeV so that this bound applies to milli-charged particles with mass 
below about 5 MeV. The upper bound on e will also apply to the model with a second 
U{1). The trapping limit e > 10~^ could be slightly modified by the presence of the 
paraphoton, but this area of parameter space is already ruled out by nucleosynthesis, 
so we do not discuss this further. In summary, SN 1987A excludes only a very narrow 
sliver of parameter space in addition to what is excluded by other arguments (Fig. 1). 

5. Milli- Charged Neutrinos 

In the Standard Model with massless neutrinos, there are four anomaly-free C/(l) 
symmetries, corresponding to the Standard Model hypercharge Ism, and {B/S — Li}. 
B is baryon number, and Li are the lepton numbers of the three lepton families. The 
hypercharge operator can be redefined to be 



CT,weak — 



Stt 



(C^y + 5Cl), 



(4.19) 




(4.20) 




(5.1) 
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without making the theory anomalous. This gives electric charge ej to z/j and gener- 
ates a proton-electron charge difference eg. Constraints on this model and variants 
have been discussed in [15, 16, 17]. 

It is more difficult to give charge to massive neutrinos in this way, as discussed 
in [30]. The solar and atmospheric neutrino deficits can be explained by masses, 
mixing the three Standard Model neutrinos. If this is the case, lepton number and 
the lepton flavours are not conserved, so the neutrinos cannot acquire electric charge 
by the redefinition of equation (5.1). If there are three additional gauge singlet "right- 
handed" neutrinos without Maiorana masses, the neutrino masses can conserve B — L 
which would allow the three flavours to have the same charge e. In this case, the 
observed neutrality of matter [1] implies e < 10~^^ [15]. It would be possible for 
v^j^ and Vj- to have significantly larger charges than this if the observed solar and 
atmospheric neutrino deficits are not due to a neutrino mass matrix mixing the three 
Standard Model neutrinos^. Bounds on the electric charge of Vr were calculated in 
[18] without assuming any relation between the tau neutrino charge and the electric 
charge of the other two neutrinos. We do not consider this possibility. 

6. Conclusion 

We have updated the bounds from astrophysics, cosmology and laboratory exper- 
iments on fermions with electric charge ee, where e < 1. These "milli-charged" 
particles could be neutrinos or new fermions from beyond-the-Standard-Model. If 
milli-charged neutrinos have mass, and the Standard Model particle content is non- 
anomalous, the electric charge of neutrinos is constrained to be < 10~^^. The updated 
bounds on milli-charged particles from beyond-the-Standard-Model are presented in 
figure 1, for both the cases where there is, and there is not, a paraphoton. For masses 
< 5 keV, a fermion with electric charge 2 x 10"^'* < ee < 1 is ruled out. An electric 
charge e > 10~® is ruled out up to masses ~ MeV. Milli-charged particles could be 
possible for masses between an MeV and a TeV. 
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